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In this paper, we use high-temperature experiments coupled with ab-initio modeling of diffu- 
sion to study An in Bi2Te3. This combined methodology enables us to elucidate key aspects of 
the mechanisms of diffusion and the solubility behavior of the Au-Bi2Te3 system. Using electron 
microscopy and energy dispersive spectroscopy, we observe fast motion of An into Bi2Te3 and con- 
centrations above the previously reported solubility limit, sufficient to deplete many microns of An 
on the Bi2Te3 surface after 165 hours at 350°C. We calculate defect formation energies and diffusion 
barriers within DFT to develop an atomistic understanding of fluxes of Au into Bi2Te3. We identify 
an interstitial mechanism responsible for the observed rapid anisotropic diffusion and provide an 
estimate of the heat of transport, Q*, associated with temperature-gradient driven diffusion. The 
calculated diffusivity along the Te-Te double layer plane is in good agreement with experiment, while 
the much slower cross-plane diffusivity requires a non-interstitial mechanism. The low formation 
energies of substitutional defects suggest these types of defects may be active during diffusion and 
may explain the high observed concentrations. 
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INTRODUCTION 

Bismuth telluride, Bi2Te3, is a narrow bandgap semiconductor commonly used for waste heat harvesting and 
electrical cooling[T] due to its large thermoelectric figure of merit, ZT. Recent interest in this material has been driven 
by the prospect of improving ZT through nanostructuring[2]-[3], as well as by the discovery of topologically protected, 
highly mobile electronic surface states in Bi2Te3 and related compounds, such as Bi2Se3 and and Bia;Sbi_a:[51 [5]. 
Both thermoelectric and topological insulator properties involve somewhat electrical transport, so understanding the 
properties of metal contacts to Bi2Te3 could facilitate further scientific and technological advances with this material. 

Bi2Te3 possesses a layered, trigonal structure based on that of tetradymite|7|i8j. Quintuple layers of Te'-^'^^-Bi-Te'-^-'- 
Bi-Te(i) are stacked in the z-dircction, as shown in Fig. [l] The mixture of covalent bonding within the quintuple 
layers and van der Waals bonding between them creates an exceptionally varied, and variable, electrical environment 
for defects and transport of atoms and electrons. Features of both bulk and surface environments are simultaneously 
present. 

Various contact metals have been used in Bi2Te3 devices, including Pt[9 , Ag lOj, and Au[TTl[T2]. Au and Au/Ti 
layers[13j form good electrical contacts with acceptable mechanical properties. In electrical devices with Au electrodes 
there is fast diffusion from the electrode into the Bi2Te3 region[T2]. Seminal measurements of diffusivity and solubility 
of Au in Bi2Te3 were carried out by Keys and DuttonTd', who found exceptionally fast diffusion in single crystal 
specimens parallel to the Te-Te double double layer planes, Au solubility (about 10^® atoms/cm'^ at 300K, much 
lower than that of commonly used dopants in thermoelectric applications), and hardening. Although the mechanism 
responsible for the high diffusivity is unknown, the diffusion mechanisms for Cu[T5, Ag[Tni[I7] and Se [TH] in Bi2Te3 
have been investigated. In the latter case, a mechanism involving antisite defects and thermal vacancies was proposed. 
In the former studies involving the metals, electron-ion interactions were thought to be relevant, partially due to the 
large carrier concentration (~ lO-'^^/cm'^) at room temperature due to intrinsic doping. Rapid diffusion in the plane 
formed by the weakly bound Te-Te double layers was suggested for Au, Cu and Ag, but there exists considerable 
uncertainty regarding the diffusion mechanism. 

Although the measured solubility is low, the importance of doping in thermoelectric [TM^ and topological 
insulator [53] applications as well as the more general question of aging in electrically contacted thermoelectric 
devices motivate us to understand the detailed atomistic mechanism of Au incorporation and transport in Bi2Te3. 
In this paper we describe experiments to observe Au diffusion into polycrystalline n- and p-type Bi2Te3. We observe 
fluxes large enough to completely deplete the several microns thick Au thin films. 

The possibility of reaching a high impurity concentration raises the prospect of significant degradation of the electri- 
cal and mechanical properties of Au-contacted Bi2Te3 over long periods of time. Furthermore, because thermoelectric 
devices are designed to be exposed to large temperature gradients, the Soret effect (migration driven by a temperature 
gradient) may also cause large concentrations of Au to accumulate in the Bi2Te3 devices. 
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FIG. 1: The 135-atom Bi2Te3 unit cell with Te atoms as small purple spheres and Bi atoms as large red spheres. The 
substitutional doping sites are labels are capitalized: Auj,j.(i) (A), AuBi (B) and Auj,g(2) (C). The three possible interstitial 
sites are labeled in lowercase: Te'^-'-Te^^' int (a), Bi-Te'^-* int (b), Bi-Te''^' int (c). The in-plane (shorter arrow) and out-of-plane 
(longer arrow) diffusion pathways are indicated by the dashed arrows in the lower portion of the supercell. 

In this paper we also compute the formation energies of isolated defects and diffusivities in order to understand 
the observed incorporation of Au in Bi2Te3. We provide an explanation of the experimental technique to observe the 
Au diffusion motivating the modeling. Through formation energies we address the question of which types of defects 
will be present. By determining a likely microscopic diffusion mechanism, we point out the origin of the observed 
anisotropic, rapid diffusion and also estimate the heat of transport. These quantities determine how quickly and how 
much the Au from a contact will be incorporated into a Bi2Te3. Our calculations suggest a second slower stage of 
diffusion associated with Au substitution in the Bi2Te3 occurs after the initial rapid diffusion measured by Keys and 
Dutton. The effective diffusivity in the second stage will be kinetically limited by the availability of substitutional 
defect sites, and thus most difficult to observe in highly crystalline stoichiometric samples. 

II. METHOD 
Experiment 

To perform aging experiments, Au contacts are sputter deposited onto two n-type Bi2Te3 substrates. The substrates 
are sintered and lightly doped with Se to a nominal composition of Bi2Te2.5Seo.5, improving the thermoelectric 
properties when compared to undoped Bi2Te3. Circular Au contacts are patterned onto the substrates for future 
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FIG. 2: The optical micrograph in a) shows a lithographically patterned Au film on N-type Bi2Te3. The Au circles range 
in diameter from 200 pm to 600 pm in steps of 50 fim. This pattern is chosen for Cox and Strack[24] contact resistance 
measurements, which will be discussed in subsequent work. The SEM micrograph in b) shows a single Au contact as processed. 
There is cutting debris on the surface because we dice a larger substrate into smaller sections. 

studies of contact resistance using the Cox and Strack method[24 . Conventional UV lithography is used to pattern 
circles ranging in diameter from 200 to 600 /im in steps of 50 /im onto the substrate surfaces. Before depositing Au, 
the substrate surfaces are briefly etched with an argon plasma, and then Au is sputter deposited to a thickness of 1 
/im as shown in Fig. [2] Aging is performed in evacuated quartz tubes at 150°C, 250°C, and 350°C for 65 hr and 163 
hr. The Au contacts are characterized with scanning electron microscopy (JEOL 7600, JEOL 840, and PhenomWorld 
desktop SEM) and scanning transmission electron microscopy (JEOL 2010 FEG operating at 200 kV). EDS analysis 
in the TEM is performed with Digiscan scan controls (Gatan, Inc.) and a SiLi EDS detector (Oxford, Inc.). 



We use a 135-atom supercell of pure Bi2Te3 for density functional theory (DFT) calculations. Previous studies of 
native defects in Bi2Te3 suggest such a supercell can accurately model isolated defects [3S]. The VASP code |26j . 
with standard projector augmented wave (PAW) psuedopotentials, with d-electrons in the core, is used to compute 
the electronic structure. A planewave energy cutoff of 500 eV and a Monkhorst-Pack k-point mesh of 2 x 2 x 2 are 
adopted. The basis vectors and ionic positions are relaxed to determine the low energy crystal structure, whose a 
— 4.35 A and c = 29.84 A lattice constants, as well as the Te-Te layer separation (3.46 A) are within 1% of the 
experimental values (a = 4.38 A and c — 30.36 A Te-Te distance — 3.63 A OUT]). Relaxations are finished when all 
forces on all atoms are less than 6 mcV/A. In the subsequent doping studies, the dimensions of the supercell are held 
fixed and only ionic relaxations are allowed. The spin-orbit interaction is taken into account because the Bi and Te 
have heavy nuclei. Fig. [T] shows the pure Bi2Te3 crystal supercell. There are three unique crystal atom sites and two 
inequivalent stable interstitial sites. 

We consider two types of substitutional defect. In one type a Au atom replaces a single Bi or Te atom and the 
replaced atom is removed entirely. In the second type a Au atom occupies a Bi or Te site and pushes the Bi or 
Te atom into an adjacent interstitial site. Interstitial defects are also modeled by adding an Au atom into the pure 
bismuth telluride supercell at one of the three possible interstitial sites with no initial displacement of any atoms. 

The method of Northrup and Zhang [55] is used to compute formation energies. For an interstitial defect 



where Eloped is the total energy of the supercell with the dopant, is the chemical potential for Au, computed 
from fee bulk Au, and Ep„re is the total energy of the Bi2Te3 supercell with no dopant. AEy is the formation energy, 
the energy penalty to create the doped system from the pure crystals. 



Modeling 



AEf — Edoped — {l^-Au + Epure) 



(1) 
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In the case of a substitutional defect, we assume the displaced host crystal atom (Bi or Te) is sent to a chemically 
well-defined reservoir. Then, 

AEf = E doped — i^J-Au + Epure — l^'Te/Bi) (2) 

where Hxe/Bi is the chemical potential of the displaced Bi or Te atom in the reservoir. If the displaced atom remains 
at the adjacent interstitial site we revert to Eq. 1. 

In a typical Bi2Te3 application, such as in the active regions of a thermocouple, there could be a variety of chemical 
environments and hence chemical potentials. Because the Au-Bi-Te phase diagram is known 23120] i we assume the 
system of Bi2Te3 with the Au dopant is in equilibrium with pure Au (from a contact, say) and either AuTe2 or BiTe, 
depending on whether the environment is Te-rich or Bi-rich, respectively. 

AuTe2 crystalizes in the calaverite structure with a three atom unit cell[3T]. The BiTe[32] structure is related to 
the tetradymite Bi2Te3 structure, with an additional Bi double layer inserted between every other Te'^^^-Te'^^-' double 
layer. There are many other Bi-rich bismuth tellurides [33]; we chose to use BiTe because both AuTe2 and BiTe 
coexist with Bi2Te3 in region of the Au-Bi-Te phase diagram between pure Bi2Tc3 to pure Au. 

The three hau, fJ-Bi, ^J■Te can be determined from the total energies of Au, AuTc2, BiTe, and Bi2Te3. We solve the 
first along with either the second or the third of the equations below. 

2^Bi + S^Te — f^BiiTes (3) 
fJ'Bi + UTe = fJ-BiTe (4) 
^^Au + 2MTe — fJ-AuTe2 (5) 

where fiAu, I^AuTe^ ^J'BiTe and ^,Bi2Te3 are the total energies of one formula unit of Au, AuTe2, BiTe, or Bi2Te3, 
respectively. In the Bi-rich case, Eq. 3 and Eq. 4 are solved, while in the Te-rich case Eq. 3 and Eq. 5 are solved. 

Since there is no ternary phase, nor a relevant BiAu phase, our choice of chemical reservoirs is nearly unique for 
the dilutely Au doped stoichiometric Bi2Te3. We note above 116 °C, BiAu2 is entropically stabilized and so it 
may appear as a possible reservoir. Errors associated with neglecting this phase would likely be much smaller in 
real applications than those due to departures from pure Bi2Te3 stoichiometry, ternary dopants for altering the host 
electronic structure, and non-equilibrium conditions. 

The activation energy for diffusion of Au in the in-plane direction between adjacent Te-Te double layer interstitial 
sites is calculated with the nudged elastic band method (NEB)[3S]. Au can diffuse by other ways in bismuth telluride. 
Along with diffusion in the cross-plane or z-direction, more complicated mechanisms, such as vacancy mediated 
diffusion, multiple-defect diffusion or even diffusion along grain boundaries in polycrystalline samples may contribute 
to the experimentally measured diffusion. The calculations presented here provide a theoretical upper bound on the 
diffusion through single crystal Bi2Te3 along solely interstitial paths. 

The starting point for the NEB calculations is an Au atom in the Te-Te double layer interstitial site in the 135-atom 
supercell. The final configuration has the Au atom moved to an adjacent Te-Te double layer interstitial site, either in 
the same plane or in an adjacent Te-Te double layer. Initially, a series of 16 equally-spaced atomic positions for the 
Au along a line connecting the initial and final Au positions are the assumed diffusion pathway. Then, forces on each 
of the 16 images of the isolated Au atom along the line are computed. All atoms in each of the 16 supercells relax, 
with the caveat that an additional fictitious force acts on each Au atom. This elastic band force is due to adjacent 
images of the Au atoms, keeping them separated. In this way, the Au image atoms are forbidden from simply relaxing 
down to the globally favored initial or final configuration sites. The forces on all the atoms, including the fictitious 
elastic band forces, in each of the image supercells are recomputed until self-consistent cycle is complete. In this way 
low-energy metastable saddle point-crossing pathways can be found in the vicinity of the initial guess for the diffusion 
pathway. We note that pathways involving substitutional defects will not be discovered using this method. 

III. RESULTS 

We observe large fiuxes of Au into the Bi2Te3 substrate, sufficient to deplete the contact pads (see Fig. [3]). The 
measured concentration is several times larger than the previously reported solubility limit [14J and even further above 
the intermediate-time net concentrations computed using the reported diffusivity (see Fig 111). The concentration of 
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FIG. 3: The figure shows a series of Au contacts after various aging conditions (time and temperature). The scale bar in each 
image is 50 fim long. The Au contacts are visually intact at 150°C up to 163 hours. At 250°C there is evidence of degradation 
around the edges of the sample after 65 hours and roughening and degradation throughout by 163 hours. At 350°C the mostly 
gone by 63 hours and is practically gone after 163 hours.. Small islands remain on the surface, and EDS suggests these islands 
are not just pure Au, but contain Bi and Te also. At intermediate aging conditions the edges of the contact seem to disappear 
first. The contact develops a swiss cheese appearance with holes that open up in the middle of the contact. This is most easily 
visible in the image after 65 hr at 350°C. 



Au in the substrate is computed by dividing the total volume of gold by the substrate volume. We fabricate two 
samples each with many contact pads, for total Au volumes of 4.18 x 10~^cm'^ and 3.55 x 10~^cm'^, corresponding 
to 2.47 X 10^^ and 5.12 x 10^* Au atoms. The substrates are both disks 1.6 cm in diameter and 0.2 cm thick. Thus, 
the Au concentration in the substrate after the Au contact pads are entirely depleted is 6.13 x 10^* atoms/cm'^and 
5.21 X 10^^ atoms/cm"^. The greater value sets a lower limit on the solubility of Au in n-type Bi2Te3 higher than that 
measured in |14j . 

The STEM and EDS analysis of the contact/substrate interface shown in Fig. [5] shows changes as the Au contact 
diffuses into the substrate. Cavities form in both the Au and Bi2Te3, suggesting both substrate and contact atoms 
migrate. The SEM images in Fig. |3] show pits (dark regions) in the Bi2Te3 substrate both under the contact and 
away from the contact at high temperatures and large times. This suggests the departure of Te or Bi atoms from 
the substrate may not be directly due to the Au contact; the Te vapor pressure of Bi2Te3 is exceptionally high and 
sublimation of Te at elevated temperatures has previously been observed. However, the EDS analysis of precipitates 
on top of the contact shows both Te and Bi migrate through the contact and precipitate. This suggests some reaction 
with the Au occurs. The calculations discussed next suggest substitutional Au defects are energetically cheap and 
the formation of substitutional Au should be kinetically limited by the cost to displace substrate Bi or Te atoms into 
interstitial sites. In the present Se-doped polycrystalline sample, there may be many Bi or Te vacancies or other 
defects whose energy can be lowered by incorporating Au into the Bi2Te3. 

Table 1 gives defect formation energies for interstitial, substitutional, and substitutional with displaced Bi or Te 
atoms. The total energy for the 135-atom pure Bi2Te3 model is -529.92 eV and the total energy for the bulk Au 
1-atom unit cell is -3.19 eV. The lowest interstitial formation energy is the Te-Te double layer site, Te^^-'-Te^^-' int, 
because of the relatively open structure and weak bonding between the Te-Te double layers compared to the other 
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FIG. 4: The parameters provided in|14) are applied to the present experimental set up to predict the time evolution of the 



concentration profile, C{x,t) = Coerfc(- 



The Au flux predicted with these parameters is much less than we observed. 



possible sites. Another interstitial site, the Bi-Te*^"'^^ int site, is unstable for Au doping. At this site the Au atom 
relaxes into the region between Te-Te double layer to the Te^^^-Te*^^^ int site. The Bi-Te*^^^ interstitial site formation 
energy is higher because the Au atom must distort more of the rigid covalent bonds present in the quintuple layer. 

The simple substitutional formation energies are lower than either interstitial defect formation energy, indicating 
these should be more prevalent. The calculated values depend on the Bi and Te chemical potentials, but are small for 
both Bi-rich and Te-rich conditions. The exceptionally low formation energy for these substitutional defects might 
be expected after inspecting the Au-Bi-Te ternary phase diagram 301, in which several Bi-rich phases coexist with 
AuTe2 and bulk Au, hinting that in the presence of Au, the Te is nearly unstable in the stoichiometric Bi2Te3. 
Recent calculations on the phonon bandstructure of AuBiTe2 36, suggest this compound is only weakly stable and 
experiments [37] show Te concentration in near-stoichiometric Bi2Te3 materials can be easily controlled through Te 
vapor pressure due to a low Te sublimation energy. This poses a problem when Bi2Te3 is annealed at high temperatures 
and the Te sublimes. 

Considering the formation energies for the substitutional models with an internally displaced Te atom, their large 
magnitude can be understood because the displaced Te atom relaxes into the Te-Te double layer disrupting the van der 
Waals bonding between the Te-Te layers. The configuration with the Au substituting for a Bi and displacing the Bi 
into an adjacent interstitial site docs not exist because it is unstable; the Bi atom returns to its crystal site and the Au 
relaxes into the interstitial region. The large formation energies of these substitutional-interstitial defect complexes 
suggests the occupation of the substitutional sites may be kinetically limited. In other words, the formation of an 
interstitial defect will be necessary in order to have an Au atom at a substitutional site, so even though substitution 
is energetically favored overall, it may take a very long time for the Au to reach the substitutional sites because of 
the large energy barrier. This kinetic mechanism also favors the interstitial pathway for diffusion, especially in highly 
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FIG. 5: In a) a HAADF STEM image shows a cross-sectional view of An on Bi2Te3. Two layers of Pt are visible atop the 
1 /im thick polycrystalline Au contact pad. The inset image is a higher magnification view of the boxed region. In it, there 
is a feature with different contrast on the surface of the Au at a grain boundary junction. EDS analysis on this region and 
the neighboring Au shows that this feature is a second-phase precipitate crust containing Bi, Te, Se, and Au. It is likely that 
atoms from the substrate diffused along the Au grain boundaries and formed a precipitate at the surface. In b), a second phase 
precipitate is seen in the middle of the Au film. Again, EDS analysis shows this precipitate contains Au, Bi, Te, and Se. There 
are also large pores (black) in the substrate and contact. 

crystalline samples, because while the formation energy for interstitials is large compared to substitutional defects, 
the effective formation energy for substitutional sites will be larger. 

The calculations yield the formation energies for Au dopants in the perfectly crystalline Bi2Te3 environment, but 
Au-dopants will coexist alongside native defects in real Bi2Te3. Antisite defects are expected to be dominant, but 
vacancies and self-interstitials may be present [23 ESj • These various native defects will allow for different Au doping 
mechanisms because of the relatively small Au interstitial and substitutional defect formation energies. The processing 
involving deformation and intentional p- or n-type doping [5^ could further alter the population and solubility of Au 
dopants. Diverse Cu doping locations in Bi2Se3, a closely related material with an identical crystal structure and 
a van der Waals bound Se-Se double layer [40], demonstrate a sensitive dependence on processing conditions in this 
system, too. Because the diffusion of Au from a contact is inherently non-equilibrium, the diffusing Au atoms in any 
case should not be expect to exclusively occupy the lowest energy substitutional sites. 

In spite of the potential for many coexisting Au doping sites, the weak bonding and large spacing between the 
Te-Te double layers suggest Au diffusion through Te-Te interstitial sites in Bi2Te3 could be important over a range of 
doping and processing conditions. Experimental measurements of diffusion barriers of Au |14j, Ag and Cu [15] 
at room temperature in single crystal Bi2Te3 indicate the diffusion is strongly anisotropic, with an activation energy 
of 0.20 eV (0.45, 0.21 eV) for the in-plane direction and 0.56 eV (0.92, 0.80 eV) in the cross plane direction for Au 
(Ag,Cu). 

We compute the Au in-plane diffusion barrier and find it is 0.18 eV, in good agreement with the experimental 
value of 0.20 eV [14]. The diffusion is not along a straight line connecting the adjacent interstitial sites. In the 
relatively open structure of the Te-Te double layer, the Au atom moves along a curved path to maximize the Te-Au 
distance. Fig [6] shows the pathway, viewed along the c-axis with only the equilibrium positions of the Te atoms in 
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Site 


Formation Energy (eV) 




Te(i)-Te(i) int 


1.0 eV 




Bi-Te(2' int 


1.13 eV 


Te-rich condition 




0.46 




AuT^eCl) 


0.08 




AU'r„(2'i 


-0.09 


Bi-rich condition 




0.14 






0.28 




Au7.g(2) 


0.11 




AuTg(i)-Te int 


3.00 




AuTg(2)-Te int 


3.42 



TABLE I: Single Au-dopant formation energies. Substitutional formation energies are computed either by displacing substituted 
lattice atom into an adjacent interstitial site (or Auj,^(x)) or by using jj,Bi and fj,Te derived assuming the system contacts 
reservoirs of Au and either AuTe2 or BiTe, in the Te rich and Bi rich cases, respectively. 




FIG. 6: Schematic illustrating the Au diffusion pathway in the Te-Te double layer plane. The Te atoms in the upper plane are 
shown as blue (color online) spheres and those in the lower plane are shown with red crosses. The Au atom is shown at the 
initial and final positions in lime green. The pathway bends upward out of the plane at its midpoint to avoid the Te atom in 
the lower plane. 
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FIG. 7: Total energy as a function of position along the diffusion pathway. The initial Te-Te interstitial site corresponds to 
reaction coordinate 0. The energies are not symmetric about the midpoint because of the underlying asymmetry in the crystal. 
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the relevant double layer shown. The Te atoms slightly distort to allow the Au atom to squeeze past, with the largest 
distortions about one-quarter and about three-quarters of the way along the transit. Fig. [7] shows the total energy 
change along the diffusion pathway. The good agreement between experimental and theoretical values of the in-plane 
activation energy and diffusivity (in the next section) suggest the in-plane diffusion in highly crystalline samples is 
indeed through the interstitial mechanism discussed here. 

We also compute the activation energy for diffusion in the out-of-plane direction to reach an adjacent Te-Tc double 
layer. In this case the barrier is nearly f .5 eV, occurring where the Au atom leaves the Te-Te double layer and pushes 
into the interior of the quintuple layer. The large barrier and the associated disagreement with the experimentally 
measured value of 0.56 eV[T3j suggest other mechanisms must allow the Au atom to move in the out-of-plane direction 
more easily. Given the low formation energies for substitutional defects in both Bi- and Te-rich conditions, diffusion 
pathways involving substitutional Au intermediate states could be involved, either through vacancies or antisite 
defects. If such pathways are relevant they would be expected to lower the diffusion barrier because the substitutional 
defect would decrease the strain in the lattice. The fact that the calculated energy barrier in the cross plane direction 
occurs as the Au is leaving the Te-Te double layer suggests other interstitial-type paths between Te-Te double layers 
will encounter a similar barrier, so some pathway not accessible to the interstitial NEB method seems to be required. 

In addition to the barrier height, the diffusivity itself D = Dqc^'"'*^'''^'^ for the Au diffusing through the Bi2Te3 
along the Te-Te double layer can be computed. Zener[41j begins from the random walk problem and suggests that 
Da = Xal^e^^/'"'. Here A is a geometrical factor taking into account the various paths the diffusing particle can 
take, ao is the lattice constant and v is the vibrational frequency of the diffusing particle in the crystal potential field. 
In a harmonic approximation v is approximately equal to ^ E^ctl'^'tna^ where E^ct is the activation energy (barrier 
height), TO is the Au atom mass and a is the jump distance. AS* is an entropy term to account for the changing strain 
in the lattice due to the diffusing particle. It is reasonable to assume AS* is zero, given the open Bi2Te3 structure 
and the small distortion of the lattice [13 [H]. For forward diffusion in the Bi2Te3 Te-Te double layer, A = 3/2; using 
ao = 4.38 Aand v = 0.49x10^^ s"^ Then the calculated Do = 1.41x10"^ qI^siUb cm^/s. The experimental value 
is I.131q'3^x10^'^ cm^/s [14], in marked agreement with the calculated value for the in-plane interstitial mechanism, 
suggesting this mechanism dominates transport parallel to the Te-Te double layer directions (and also that AS" is 
indeed nearly zero). 

Because thermoelectric devices are usually exposed to significant temperature gradients, it is useful to consider 
diffusion under such conditions. A temperature gradient introduces an additional term in the expression for diffusive 
flux J in a concentration field, c. 

where Q* is the heat of transport. The second term on the right captures the Soret effect. We estimate Q* as 
follows. The simple interstitial transport mechanism has a relatively low activation barrier and so we assume this 
mechanism dominates transport in the double layer. Because the AS' in the diffusivity is small, and no covalent bonds 
are formed or broken, we further simplify the estimate of Q* by assuming the phonon spectrum along the diffusion 
path is unchanging. In this approximation -Eact < Q* < Eact |42] ; the limits are satisfied when the activation energy 
is dissipated about the initial site, or when all of it is instead carried to the final site. If the activation energy is 
evenly dissipated between the initial and final sites then there is no energy flux as the particle moves and Q* = 0. In 
the present case the diffusion barrier profile, with the largest drop in energy occurring near final site, suggests energy 
will be mostly dissipated near the final site and so we estimate Q* ^ 0.18 eV (the calculated diffusion barrier) for 
in-plane thermomigration. Experimental measurement of 0.265 eV heat of transport for Ag in Bi2Te3[43] suggests 
the sign and magnitude of this estimate are reasonable. 

To estimate the relative importance of the concentration- and temperature-driven (Soret) diffusion fluxes, we 
consider a model system with midpoint temperature of 200°C, a 100°/cm temperature gradient and an initially 
uniform concentration at the solubility limit . Using the solubility limit measured in |14| of 1x10^^ atoms/cm'^ and 
a diffusivity of 1.36x10"^ atoms/cm^-s the flux from the Soret term at the midpoint will be 6.37x10"'^'^ atoms/cm^-s. 
Using the same diffusivity and temperature concentration gradient-driven flux term at 1.66 hours from Fig. |4] is 
2.5x10"'^^ atoms/cm^-s, already an order of magnitude lower than the Soret term. Using instead our higher measured 
solubility limit of 6x10^^ atoms/cm^ will simply multiply both flux terms by six. This estimate suggests Au fluxes 
from contacts into Bi2Te3 operating in thermoelectric modules will be largely temperature gradient-driven because 
the initial concentration-driven diffusion occurs so quickly. 
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IV. CONCLUSION 

We observed diffusion of Au into n-typc Bi2Tc3 substrates and migration of Bi and Te on to the top of the Au 
contact material. We have calculated formation energies, diffusion barriers and diffusivity for a likely path for single 
Au atoms diffusing in Bi2Te3. The lowest energy interstitial defect is between the Te-Te double layer. The other 
stable interstitial site is in the quintuple layer between the Bi and Te layer. Substitutional defects have low formation 
energies for Bi2Te3 in equilibrium with Au and either AuTe2 or BiTe chemical reservoirs. Au substitutional defects 
with an adjacent displaced Bi or Te atom are found to have larger formation energies than interstitial defects due to 
the presence of, essentially, two defects. 

Calculated values for the diffusion barrier and diffusivity along the in-plane direction agree well with experiment. 
The diffusion pathway winds through the Te atoms in the double layer. For diffusion in the perpendicular direction, 
the barrier is large and the agreement between the calculated value and experiment is poorer. The discrepancy may 
be due to diffusion pathways beyond the reach of the NEB method we employ. 

Our results suggest that in Bi2Te3 near stoichiometry, if and when equilibrium is reached, the majority of Au atoms 
should occupy substitutional sites and any Au interstitial defects should predominately be found in the Te-Te double 
layer. The initial diffusion of Au should occur mainly along the in-plane direction by a mechanism involving the 
intertstitial sites between the Te-Te double layer. When a diffusing Au atom occupies a low-energy substitutional 
site, it will cease moving. Because of the kinetic barrier associated with moving the displaced Te or Bi atom into 
the interstitial, however, it may be relatively rare for an Au atom to occupy a substitutional site. If an Au atom 
occupies a Te-Te double layer interstitial site, it may be expected to rapidly diffuse along the in-plane direction. 
We suggest the diffusion occurs in two stages. First, Au rapidly diffuses along the Te-Te double layers then quickly 
ceases to diffuse as the low solubility limit for interstitial Au is reached. Next, slower diffusion through an unknown 
vacancy- or interstitial-mediated process occurs, filling the substitutional sites. The timing of the second stage might 
be expected to depend sensitively on the composition of the Bi2Te3; for high quality single crystal specimens it would 
be exceedingly slow because of a lack of vacancies and other defects to aid the formation of substitutional Au defects, 
while for more highly defective poly crystalline materials the second stage would occur more quickly. 

Our experiment on Au in contact with Bi2Te3 demonstrates the rapid diffusion of Au continuing beyond the 
solubility limit suggested in [T3]. From a practical perspective, some of the negative consequences of Au diffusion 
from contacts into a thermoelectric device could be ameliorated by increasing the volume of the Au contacts so they will 
not be completely depleted. However, voiding may still occur at the interface (as we observe), substantially reducing 
the effective contact area. More concerning is that the incorporation of Au into the active thermoelectric material 
may gradually spoil the carefully tuned doping and transport properties. While for highly crystalline materials the 
second stage substitutional diffusion is expected to be kinctically slowed by the large cost of forming host interstitial - 
Au substitutional defect complexes, these materials might also be expected to show larger deviations in thermoelectric 
performance due to incorporation of Au than would more polycrystalline or off-stoichiometric samples. 
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